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Abstract: We use AdS/CFT to investigate i) high energy cohisions with balls of de- 
confined plasma surrounded by a confining phase and ii) the rapid localized heating of a 
deconfined plasma. Both of these processes are dual to collisions with black holes, where 
they result in the nucleation of a new "arm" of the horizon reaching out in the direction of 
the incident object. We study the resulting non-equilibrium dynamics in a universal limit 
of the gravitational physics which may indicate universal behavior of deconfined plasmas 
at large Nc- Process (z) produces "virtual" arms of the plasma ball, while process (ii) can 
nucleate surprisingly large bubbles of a higher temperature phase. 
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1. Introduction 

Studies of gravity/gauge duality |jl], ||, ^ have given insight into many aspects of strongly 
coupled gauge theories. A central aspect of this correspondence is that deconfined plasmas 
in the gauge theory are dual to black holes on the gravity side of the correspondence Q. 
Although no gravity dual is known for gauge theories of direct experimental interest (e.g., 
QCD), the universal properties of black holes in gravity suggest that insights obtained from 
gravity/gauge duality may describe universal aspects of deconfined plasmas at large A, Nc, 
where Nc is the number of colors and A is the 't Hooft coupling. The known utility of the 
large Nc expansion in QCD then further suggests that such dualities may help to explain 
the physics observed at heavy ion colliders 0, |5[ |6|, 0. See e.g., |8| for a recent review. 

Much of the literature on this subject has focused on using black holes to explore equi- 
librium thermodynamic properties of deconfined plasmas (such as entropies, free energies, 
and phase transitions, see e.g., ||3|, |9|) or quasi-static near-equilibrium properties such as 
hydrodynamic transport coefficients, see e.g., |5|, While transport coefficients can de- 
scribe important hydrodynamics (such as elliptic flow, see e.g., fl^), such processes must 
be quasi-static in the sense that each local region of the plasma remains close to thermal 
equilibrium. 

However, one may also ask about more dynamic processes. One example is the forma- 
tion of a plasma from collisions studied in 0. We will be interested here in related (but 
different) settings, involving the rapid transfer of energy to pre-existing plasmas. At low 
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rates of energy transfer, these processes result in a flow of heat described by hydrodynamic 
thermal conductivity, perhaps associated with some expansion or hydrodynamic flow of 
the plasma. But what happens at higher rates of energy transfer? 

The result is most easily seen from the gravity dual. Recall that the most interesting 
feature of the dual black hole is its event horizon and that this horizon is a surface defined 
by null geodesies. A small energy flux across the black hole horizon leads only to local 
expansion of the horizon in a manner readily described by hydrodynamics (see e.g., [pi]]). 
On the other hand, larger fluxes lead to more dramatic results. As an extreme case, consider 
the collision of a small black hole with a larger black hole as shown in Fig. ||. 

At first, the two black holes approach each other, each 
maintaining its own event horizon with only minor effects from 
the other black hole. However, when the black holes are suffi- 
ciently close together a new common horizon suddenly forms 
around both black holes, resulting in a single highly deformed 
black hole. In particular, the region between the two black 
holes is now enclosed in a new "arm" extending from the large 
black hole. The arm is then absorbed by the black hole as the 
system relaxes to equilibrium. 

This picture is currently being investigated in numerical 




simulations of black hole collisions |12], though the analogous 
phenomenon for collisions of equal mass black holes is already 
well illustrated by numerical simulations, see e.g., |l^. In any 
case, it has long been known on general grounds (see e.g., p^ , 
|T^) that the new arm of the horizon forms along a spacelike 
caustic; i.e., the new arm of the black hole effectively nucleates 
"from scratch." Furthermore, the incident object need not be 
another black hole for this to occur; any sufficiently rapid flux 
of incoming energy results in the formation of a similar black 
hole arm that reaches out toward the incoming object. As 




fl.l) 



Figure 1: This cartoon de- 
picts a time sequence of the 
creation of an "arm" as a 
small object falls into a large 
black hole. Note that the 
arm forms before the arrival 
of the incident object, 
discussed in ||16[ (which generalizes the arguments of 111]), in 

any spacetime dimension d the threshold for such nucleation to occur is^ 

E T 

where E, A are the energy and transverse area of the incident object, T is the temperature 
of the target black hole, Gd is Newton's gravitational constant, and we have set h= 1. 

The purpose of this paper is to calculate the shape of the nucleating arm and to 
interpret the results in terms of deconfined plasmas; we shall be less concerned with the 
final relaxation to equilibrium. As we discuss in more detail below, the formation of a 
black hole arm can correspond to either of two processes for the plasma, depending on 
the direction in which the arm extends. If the black hole arm reaches out along one of 
the gauge theory directions, then our process is naturally interpreted as the formation of 
a "virtual" arm of deconfined plasma catalyzed by the incident fiux of energy. If, on the 
^This result holds for a homogeneous object much smaller than the curvature scale at the black hole 
horizon whose cross-section transverse to the motion is round. The size of the object along the direction of 
motion does not affect the result so long as the collision takes place in less than the relaxation time 1/T of 
the black hole. 
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other hand, the black hole arm extends in the holographic direction (which corresponds to 
energy scales in the gauge theory) , then it represents a strong local heating of the plasma 
via, for example, the decay of an excited quasi-particle of large mass^. Such heating takes 
the plasma far from local equilibrium and is not well-described by hydrodynamics. In 
particular, when the system allows further phase transitions, local heating can create a 
bubble of a higher temperature phase within the original plasma. 

Though the details of specific models are certain to contain much interesting physics, 
our goal here is to take only the first steps toward understanding this dynamics in the 
context of gauge/gravity duality. For this reason, we focus here on universal or near- 
universal properties of collisions with black holes (and thus presumably with deconfined 
plasmas). In particular, we use the fact that the geometry sufficiently close to the horizon 
of any black hole of non-zero temperature can be approximated by Rindler space, which 
is just flat Minkowski space in accelerated coordinates. We also take the incident object 
to be both small and highly boosted so that horizon dynamics can be studied through the 
simple exercise of tracing null geodesies through the Aichelburg-Sexl metric |17]. 



This analysis is performed in section y. We then interpret the results in terms of large 
Nc gauge theories in sections ^ and |^. Section |3| considers high energy collisions with plasma 
balls and the resulting virtual arms that extend in the gauge theory directions. Section 
^ considers the case of rapid local heating. As a particular example we examine certain 
gauge theories with flavor in which stable mesons can exist within the original plasma. In 
that case, our local heating can cause the mesons to melt within a localized bubble, though 
the bubbles turn out to be surprisingly large. Gross properties of such bubbles are studied 
at both strong and weak 't Hooft couplings. We conclude with a brief discussion in section 

i 

2. How far can a black hole reach? 

In this section we address universal properties of collisions with black holes. Our primary 
goal is to characterize the black hole arm that nucleates to engulf an incident object. 
We wish to determine the shape of this arm, and in particular its thickness and length. 



Section 2.1 discusses the physical setup and determines the resulting horizon, while section 



2.2| extracts the desired parameters. 



2.1 Physical and mathematical framework 

Our study of collisions with black holes is carried out under two key approximations, 
within which the behavior is universal. First, we use the fact that the geometry sufficiently 
close to the horizon of any black hole with temperature T ^ can be approximated by 
Rindler space, which is just flat Minkowski space in accelerated coordinates. Second, we 
consider the limit where the object approaches the black hole at high velocity so that generic 
(uncharged, non-spinning) objects are well-described by the simple Aichelburg-Sexl metric 



^One can see from the gravity dual that such quasi-particles exist at large Nc,\ even in the decon- 
fined phased. These quasi-particles correspond to bulk gravitons, strings, or localized black holes raised 
sufficiently far above the horizon of the black hole associated with the deconfinement transition. 
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. Studying horizon dynamics during the colhsion then reduces to tracking nuh geodesies 
(the so-cahed generators of the event horizon) as they propagate through the associated 
spacetime. Such geodesies were studied in Q|, which we review below for completeness. 




Figure 2: Left: Null geodesies (thiek black lines) are shown projected into the Z, X° plane. The 
geodesies experience a discontinuity when they reach the plane 1^ = Vo of the Aichelburg-Sexl 
wave, shown as a thin line. The dashed lines are the U — and ^ = planes. In this example, 
the velocity component dZ/dX'^ of the geodesic is positive both before and after the shift, but in 
general for y < Vb the velocity dZ/dX^ may be positive, negative, or zero depending on b. Right: 
Null geodesies are shown projected into the plane. Below V ~ Vq, the geodesies focus 

toward X^ = 0. Both figures show b = 1.2, c^j/i =1, Vq — 2, and d = 5. 



We consider the general case of d > 4 spacetime dimensions. In light-cone coordinates 
V = + Z, U = X^ — Z the Aichelburg-Sexl metric |l7[ takes the form 



ds^ = -dUdV + dp' + p^dJ7^_3 + ^iHpW " VQ)dV' , (2.1) 

where we have defined p'^ = {X^)^ + . . . + (X'^"^)^. Here the incident mass moves in the 
negative Z direction at p = within the null plane V = Vq. Note that, due to the delta- 
function above, the metric differs from flat Minkowski space only within this null plane. 



The result (2.1) can be obtained from the mass m Schwarzschild metric in d dimensions 
by Lorentz transforming to a frame in which the particle travels with speed v. In units 
where c = 1, one then takes the limits m ^ and v ^ 1 in such a way that the quantity 
p, = m(l — t;^)~^/^ appearing in ( p.ip remains constant. 

The Aichelburg-Sexl potential'^ appearing in ( |2.l[) is given by 



Hp) 



-C4 In if d = 4 

2.2 

^ ifd>4 ^ ^ 



(d-4)p' 



^In contrast to the usual convention, we have separated the factor of /i from the rest of the potential in 
order to make "l'(p) boost-invariant. 
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= 7^ , (2.3) 

On denotes the volume of 5""", and the length scale po is an arbitrary gauge choice. Note 
that the potential ^*(/o) solves Poisson's equation in d — 2 dimensions with a point source. 

The Aichelburg-Sexl solution will accurately describe the metric for p'> r, where r is 
the transverse size of the incident object. The metric for p < r depends on the internal 
structure of the object; we will ignore the details in this region and treat p ~ r only as a 
cut-off on the singularity of <I>. To describe a beam of incident objects, this point source 
would be replaced by a line, or by a line segment for a pulsed beam. We consider here only 
the point source case (2.2), or equivalently the limit of short pulse duration At <C 1/T. 



In the Rindler approximation, the black hole time-translation is the boost generator 
X = 2-KT{Vdv — Udu), where we have used the fact that the so-called surface gravity 
K which determines the normalization of the boost generator is related to the black hole 
temperature through k = 2irT in units with h = c = 1. Thus the energy of the incident 
object is 

E = -XuP" = 2TipTVo. (2.4) 
We wish to use the above metric to locate the (future) black hole event horizon, which 



as always is determined by a boundary condition in the far future. Since (2.1) is precisely 
the Minkowski metric for V > Vq^ the horizon there must agree with some null plane; i.e., 
with a Rindler horizon. We take this null plane to be C/ = 0. The past horizon of the 
black hole is similarly determined by a boundary condition in the far past; we take this 
horizon to lie at ^ = 0. Having established the boundary conditions, we need only trace 
the relevant null geodesies back from the far future through the shock wave at y = Vq to 
find the full horizon. We may think of this as calculating the deviation of the true event 
horizon from the unperturbed event horizon at U = 0. 
The geodesic equation 

was studied for (|2.1|) in |jl^ and it is instructive to review their calculation. Setting the 



index v = V ui (2.5) reveals that V is an affine parameter, so henceforth we will set \ = V. 

(Pp p, 



Setting V = p then yields 



dV^ 2 dp 



5{V - Vo) = 0, (2.6) 



while V = U leads to 

g,-,M,)±m-V,)-2/^^piV-V„)=>>. (2,7) 

These equations are solved by 

p{V) = b + p^{Vo-V)Q{VQ-V), (2.8) 
U{V) = p^^T^iV - Vo)Q{Vq -V)- p^{b)Q{Vo - V) , (2.9) 
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where we have chosen boundary conditions such that far in the future (i.e., V > Vq) we 
have p — > 6, [/ — > (see Fig. |2|). It is straightforward to check that ds^ = along such 
curves, so these are indeed nuh geodesies. 

For V > Vq, the geodesies are straight lines of x° 
constant p in the U = plane. When a geodesic 
of given impact parameter p = b reaches V = Vq, 
there is a discontinuous jump^ in U given by AU = 
— ;U$(6), where the sign corresponds to tracing the 
geodesies backwards from V = oo. For d > 4, we 
have AU < 0, and AU as b ^ oo. This is not 
the case for d = 4 due to a logarithmic divergence. 

Below the discontinuity our null geodesies are 
once again straight lines, but since 
dp $'(6) 



1 + 



> for F < Vo, (2.10) 




they now begin to focus toward p = (Fig. ^). 
All geodesies with impact parameter b intersect at 
a common point, and the set of such intersections 
traces out a caustic as a function of b. Physically, 
this caustic describes the nucleation of the event 
horizon. 

From (|2.8| ) and ( |2.9| ) we see that the caustic oc- 
curs at 

Vc = Vo + 



Figure 3: The caustic is shown in the 
Z,X^ plane for Cdfi = I, Vq = 2, and 
c? = 5. For b > b^:, the caustic is hidden 
behind the past horizon. For b < b^, 
(black segment) the caustic lies above 
the past horizon. This portion of the 
caustic nucleates the "arm" of the black 
hole. 



26 



p^'{b) 



U, 



-/.$(6) + ^$'(6) 



(2.11) 



(2.12) 



This curve is plotted in Fig. ^. We have distinguished the F < and F > parts of 
the caustic, as only the V > Q part of the spacetime is physically relevant; the rest is 
hidden behind the past horizon. Since geodesies with large b emerge from the past horizon, 
we will refer to them as forming the "body" of the event horizon, while geodesies with 
small b intersect the caustic at y > and form the "arm" of the horizon described in the 
introduction. We see from Fig. ^ that the caustic along which the arm nucleates approaches 
the world line of the incident object in the limit 6 — > 0. Thus the arm reaches out in the 
direction from which the object is incident by an amount that depends on the cut-off r. 

2.2 Extraction of parameters 

Within our approximations, the above equations determine the shape of the horizon and 

"^A p-independent relative shift of the regions above and below V = Vo is a. gauge transformation which 
changes AC/ by a constant. It is in this sense that po in (2.2) is a gauge parameter for d = 4. In contrast, 
the p-dependent part of At/ is a physical effect describing the integrated relative shift of geodesies 
passing through the shock wave. Were the delta function in ( ^.l| ) replaced by a smooth function of compact 
support, the geodesies would become continuous but would still experience a net displacement A[/ as they 
traverse the shock. 
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Figure 4: Left: The distortion of the body of the horizon by the gravitational field of the incoming 
object. The curve shows where geodesies with 6 > 6, intersect the past horizon. The cusp corre- 
sponds to 6 = 6*, where the caustic intersects the past horizon. Right: The body {b > 6*, lighter 
part of the curved surface) and arm (6 < 6^ , darker part of the curved surface) of the horizon. The 
surface shown is generated by geodesies for a range of impact parameters b. The portion of the 
surface above the caustic and with < Vq is shown. For small b the spacelike curve of caustics 
from which the arm nucleates is seen to approach the line {V — Va, — 0) showing the path of 
the incoming particle. Both figures show CdfJ. — I, Vq — 2, and d = 5. 



the newly-nucleated arm. We now extract parameters describing the various length scales 
involved for use in our discussion of deconfined plasmas. Since both body and arm are 
dynamical objects, their lengths (and widths) will be functions of time or, more generally, 
of the spacetime surface on which they are measured. We will make certain natural choices 
below. 

Let us begin with the body of the future event horizon, which consists of geodesies that 
emerge from the past horizon. The point where a given geodesic intersects the unperturbed 



past horizon is given by setting y = in (2^) and (|2.9| ). The corresponding curve is plotted 



in Fig. ^ (left), and shows the distortion (at the past horizon) of the black hole body by 
the incoming flux of energy. For d > 4, geodesies with large b experience only a small 
deflection, and the body at large b approaches the unperturbed event horizon. The cusp in 
Fig. ^ (left) is the event at which the caustic emerges from the past horizon. It is described 
by K = in ( p. 11 ) and thus by impact parameter^ 



i_ 

-2 



Both the evolution of the body and the nucleation of the arm can be seen (Fig. ^ (right)) 
by plotting the geodesies for many values of b. 

We now extract some relevant length scales describing the deformation of the unper- 
turbed horizon, including the length of the arm. One useful measure of distance is the 

^Tlie result &* ~ (E/T) <i-2 is quite general and does not depend on the use of the Aichelburg-Sexl 
approximation. See p^ ]. 
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proper distance from the surface U = 0,V = where the unperturbed horizons intersect. 
Surfaces of constant proper distance contain the worldhnes of static observers outside the 
unperturbed black hole. These surfaces are just hyperbolae in the Minkowski coordinates 
{X°,Z) satisfying 

^2 _ ^^0^2 ^ ^ ^ ^2.14) 

for some constant a. 

We can ask which of these hyperbolae intersect a geodesic with given impact parameter 
6, though we must take some care due to the jump across the Aichelburg-Sexl plane. Let 
P denote the point on the geodesic immediately below the shift, so that U{P) = —fi^{b), 
V{P) = Vq. There are two cases: In the first case, dZ/dX^ at P is greater for the geodesic 
than for a surface of constant a. In this case the geodesic is at a < a{P) for V <Vq and we 
may say that this part of the horizon extends precisely to a[P). This occurs at large 5, say 
b > b. In the second case (b < b), one finds that dZ/dX^ at P is less for the geodesic than 
for a surface of constant a. In this case, the extension of the horizon in the Z direction is 
characterized by the value a associated with the hyperbola tangent to the geodesic below 
P. The two cases are shown in Fig. |5|. 




Figure 5: Null geodesies and some hyperbolae for constant a are shown for Cd — 1.5, Vb = 4, and 
d = 5. Case 1: For b sufficiently large (e.g., b = 1.3), the desired hyperbola intersects the geodesic 
at the point P, that is, immediately after the shift. Case 2: If 6 is smaller than some critical value 
(e.g., for b = 0.95), the geodesic is tangent to the desired hyperbola at a point with V <Vo. 

We restrict our analysis to the case d > A, where the perturbed horizon asymptotically 
approaches [/ = so that there is a well-defined asymptotic bifurcation surface from which 
to measure proper distance. The results are: 



Case 1 : For b > b a = a{P) 



{d-md-3b^-^T' 



Case2: For 6 < 6 « = ^ + (2-15) 
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where b is given by 



1 , , 1 



The maximum extension (or length) of the arm^ is just ( 2.15D evaluated at the cut-off 



b = r. Note that for d > 6 any geodesic intersecting the arm satisfies 6 < 5* < 5, so that 
case 2 always applies. Even for d = 5, case 2 applies whenever the effect is significant, 
r <^b^. In this limit one finds 

LbH arm ^ f^f^^- (2-17) 

We note once again that, due to the Rindler approximation, this result is justified only 
in cases where Lbh arm is smaller than the characteristic curvature scale of the spacetime 
near the horizon. 

We may also characterize the arm by its width in the Aichelburg-Sexl plane. At each 

1 

4:GdE \ d-2 



U this width is just the impact parameter, which ranges from 6^, = yoj^^j j where the 
arm connects to the body, to the cut-off r at the tip. 

3. Arm nucleation for plasma balls 

In gauge/gravity duality, balls of deconfined plasma surrounded by a confining phase are 
dual to black holes that are localized in the gauge theory directions. Such black holes 
naturally have temperature T = T^, where is the temperature of the deconfinement 



phase transition |20|. A particle incident on the gauge theory plasma ball corresponds to 
a flux of energy incident on the black hole. If this flux is large enough, the black hole may 
nucleate an arm as described in section |2|. 

How precisely should we understand the 

corresponding effect in the gauge theory? In 
some sense the plasma ball must also nucleate 

an arm along the path of the incoming parti- Figure 6: A particle incident on a black hole 
cle. What is interesting is that Fig. | (right) ^ual to a ball of deconfined plasma, 
indicates that this nucleation occurs in front of 
the incoming particle; i.e., before the particle arrives. 

Such behavior may appear to be acausal, but the same words might be applied to 
the nucleation of the arm on the black hole. Yet the physics of the black hole arm is 
clear. Gravity propagates causally, and the spacetime in the region of the arm has no 
locally measurable differences from the spacetime outside. It is the fact that one traces 
the event horizon backwards from the far future that allows the arm to form before the 
particles arrive. In other words, our designation of the arm as being "part" of the black hole 
describes the fact that, due to the imminent arrival of the incident energy, the gravitational 
interactions will become large enough that no information from the arm-like region will be 
able to escape to infinity. 

®In the regime of most interest, r ^ b*, the contribution from the black hole body is negligible in 
comparison to the length of the black hole arm, so henceforth we will refer to the "maximum extension of 
the arm" as simply the "length of the arm." 
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This suggests a straightforward interpretation of the corresponding "arm" of the plasma 
bah: Due to causahty, the arm is locahy indistinguishable from the confining vacuum. In 
particular, its energy density remains that of the vacuum and the arm contains no plasma. 
Instead, the arm marks a region of spacetime in which any additional particle has a large 
probability to be caught up in the interactions that will ensue and thus to fuse with the 
plasma ball. In effect, the arm describes a temporary and highly anisotropic enlargement 
of the capture cross-section associated with the plasma ball. To remind the reader of this 
interpretation, we will henceforth refer to such arms as "virtual" in the context of plasma 
balls. 

Now, even without the nucleation of vir- 
tual arms, the capture cross-section for par- 
ticles incident from infinity will be somewhat 
larger than the size of the plasma ball itself. 
Recall, for example, that while the event hori- 
zon for the familiar static Schwarzschild black 
hole lies at area-radius 2G4M, any free particle 
passing within the sphere of area-radius 3G4M 
is necessarily captured by the black hole. It is 
thus clear that there will be a similar effect in 
the gauge theory dual''. 

However, the arm describes a region of space 
time tied even more strongly to the plasma ball. 
Suppose that two additional (probe) beams are 
aimed so as to collide near the plasma. If the 
interaction region lies outside the virtual arm, 
then some of the debris from the interaction 
can escape to infinity. However, when an inter- 
action takes place inside the virtual arm, all of 
the particles produced in the collision will have 
large amplitudes to fuse with the plasma ball 
and no signal will reach an external detector. 
This difference is illustrated in Fig. |^ 

One may hope to one day observe the ef- 
fects of similar virtual arms in experiments with 
deconfined QCD plasmas. In this context, re- 
call that our description of black hole arms in 
section |2| is universal in the limit of i) high velocity collisions and ii) studying the region 
near the black hole. As a result, experiments that probe plasmas in the corresponding 
regime can provide a new universal test of the extent to which such plasmas admit a 
dual description by any semi-classical gravity system, without first needing to propose a 
particular dual theory. 

^Our Rindler-space near- horizon approximation is enough to predict that any low- mass particle incident 
from infinity is captured as long as it approaches within some distance of order the curvature scale of the 
spacetime near the horizon. However, the full capture cross-section may well be even larger. 




Figure 7: Top: Particles produced in a col- 
lision exterior to the virtual arm can reach 
the detector. Bottom: Particles produced 
in a collision inside the virtual arm fuse with 
the plasma ball and will not be detected far 
away. 
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To illustrate this idea, let us consider the length of the arm ( |2.17| ) , which describes the 
limit of large velocity for the incident particles. We wish to translate this result into a length 
Lpiasma arm for the arm of virtual plasma measured in gauge theory units. Since ( 2.17| ) 



contains a factor of G^, and since the relation between Gd and the gauge theory parameters 
is model-dependent, the result will no longer be completely universal. However, the scaling 
of Lpiasma arm with various parameters describing the incident particle will be universal. 
Precise predictions can of course be obtained for particular models; we give an example at 
the end of this section. 

The main task is to determine how distances (such as r, Lbh arm) in the gravity 
description are related to gauge theory distances. Here the high symmetry of our near- 
horizon limit is quite useful. Recall that we consider regions close enough to the horizon 
to regard the surface as flat. Thus we have Euclidean symmetry along the horizon as 
well as translation symmetry in time. There is also a reflection symmetry which inverts 
any constant z surface through any point on that surface. As a result, in terms of the 
coordinates p, z, t the gauge theory metric must be of the form 

ds^M = o-^(^) {-f{z)dt^ + g{z)dz^ + dp'^ + p^dnl)_^) , (3.1) 

where we have taken the gauge theory to live in D dimensions. In addition, the original 
coordinate t was normalized to the unit time-translation in the gauge theory. Thus, / = 
o""^. Since the gauge theory space-translation (say, in the transverse directions) also has 
constant norm, a must be a constant. The gauge theory metric is thus 

ds^M = -dt^ + {9{z)dz^ + dp^ + p^dnl)_3) . (3.2) 

Any such metric describes flat Minkowski space (as desired), but the particular form of 
( |3.2| ) tells us how to translate our bulk coordinates into gauge theory terms. 

Information about g{z) may now be obtained by considering the energy of strings 
stretched in both gauge theory and bulk. First recall that ( |2.17| ) was obtained in the limit 
where one is close enough to the black hole horizon to approximate the metric by the 
Rindler metric, 

ds^ = -{2TTTdfz^dt^ + dz^ + dp^ + p'^dnj_^ . (3.3) 

Now consider a bit of string located near the horizon and oriented parallel to the horizon, 
that is, at constant z. If, for example, such a string were stretched in the p direction (at 
fixed angles) with coordinate extent Ap, then its energy would be ^ = ej^TdzAp. Far 
from the plasma ball, the energy of the corresponding stretched string in the gauge theory 
would be £' = {27ra')~^{aAp) where a' is the gauge theory Regge slope. However, these 
two expressions cannot be equal as the gauge theory expression does not depend on z. The 
point here is that there is an attractive force between the string and the plasma ball, so 
that proximity to the plasma ball changes the effective string tension (and thus the effective 
Regge slope) by a ^-dependent shift in the tension. 

Let us call the effective Regge slope (3{z)a' . Setting the resulting energies equal yields 
CT _ ^ j£ assume that the renormalization factor /3 depends only on z and not on 
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the orientation of the string, we may obtain information about g{z) by considering bits of 
string stretched in the z-direction: 



£;^TdzAz = E = ^;^a{^Az), or g = I. (3.4) 

This imphes that spatial distances are merely rescaled by a factor of a when passing 
from gravity to gauge theory measurements. Thus, the size of the incident object in the 
gauge theory is f = err and the length of the virtual plasma arm is 

^^ 2G,^'^-V^ E 

^plasma arm ^ N^Q,^_^Td f'^"^ ' ^'^'"^^ 
T2n . IBd-2 ;^ - ^ ( \ 



where Gd and a are defined by Gd = N'^Gd/^i''^ and I?' = <7 I -7 I . We conclude that 
the length of the virtual arm scales with the energy of the incident particle and inversely 
with a power of the particle's characteristic size. This power determines the effective 
dimension of the gravity dual. Similarly, the width of the arm ranges from W^^^^^ = 

(^Oi'^^"^ (iv^^^Tj) "^'^ where it connects with the plasma bah to W^^^^^ = f at the 
tip. 

It is clear that Gd is a dimensionless model-dependent function; i.e., a function of A. 
To clarify the nature of a, let us first consider a = 27r (3a' Tdz/i'^. For a given plasma ball, 
we already know that o" is a dimensionless constant summarizing the scale of gauge theory 
proper distances relative to gravity proper distances near the black hole. On the other 
hand, far from the black hole the scale factor relating gauge theory relative to gravity 

distances is just [jt] ■ Thus, = a [ ] gives precisely the scale factor of gravity 



distances near the black hole relative to those far away; i.e., it is determined by the black 
hole metric. As a property of a classical solution of a pure gravity theory, a cannot depend 
on Gd, and must be a function only of L/R, where R is the size of the black hole and L 

is the AdS scale. Let us, however, consider the relation a = a = lirPa'^^'^Tdz/is 

from the gauge theory point of view. (In terms of the gauge theory proper distance z = az 
this is (T^ = 2nPTdZ.) The only factor through which the size of the black hole (i.e., the size 
of the plasma ball) might enter is /3, the renormalization of the Regge slope. But since this 
renormalization is due to interactions with the plasma ball, at small z (small separation 
from the plasma ball) one expects (3 to depend only on z and to become independent of R. 
Thus we deduce that a is in fact a constant of order one which depends only on qualitative 
features of the model. From the gravity perspective this means that a is independent of 
Gd, L, R, while from the gauge theory perspective a is independent of A, Nc, a' . 

A particular gauge/gravity duality of this kind was considered in ||3|, There a 

(2+1) dimensional theory was obtained by taking the low energy limit of = 4 super 
Yang-Mills compactified on a Scherk-Schwarz spatial circle. This theory is dual to a IIB 



supergravity system whose ground state is the Cartesian product of the AdS soliton |21| 
and S^. In the A — > limit, the field theory is pure (2+1) Yang Mills |p. At large A, where 
the classical gravity description is valid, the gauge theory is unfamiliar. Nevertheless, 
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it was conjectured in ||20[] that finite mass black holes have a dual description in this 
theory in terms of localized lumps of deconfined plasma. If the incident particle has no 
i?-symmetry charge, the corresponding object in the gravity dual is delocalized on the 
so that the effective dimension is d = 5 and (using the conventions of section |4.1| below) 
G5 = 2-V4^A3/4. 

One would also like to determine the constant a for this model. In principle, this can 
be read off from the bulk metric which describes the domain- wall limit considered in |20| 
(where the black hole fills a half-plane). In that case a = \k\oo/\k\BH , where |A:|oo is the 
norm of a translational Killing field k along the domain wall evaluated far from the black 
hole but on the IR fioor of the spacetime and |A;|B_ff is the norm of k evaluated near the 
horizon of the black hole. Using Fig. 6 of pOl, it is reasonable to estimate that ~ 1. 



4. Localized heating of deconfined plasmas 

Section ^ above discussed the gauge theory dual of a process in which a high-velocity object 
is incident on the black hole from a gauge theory direction. However, another interesting 
setting occurs when the object is incident on the black hole from the holographic direction. 
In this case, the incident object represents a localized heating of the plasma, perhaps due 
to the decay of massive quasi-particles (see section ^ below) . 

If the theory has additional phase transitions above the deconfinement temperature T^, 
then this local heating can induce the formation of a bubble of the high temperature phase. 
We now study such bubble formation and the ensuing bubble dynamics in a particular 
gauge/gravity correspondence. 



4.1 Brief review: AdS/CFT with flavor 

Recall p2| , 23] that in AdS/CFT, a small number of fiavors {Nf <^ N^,) of fundamental 



matter ("quarks") may be described by probe D-branes in a background AdS space. The 
transverse displacement of the probe branes from the D3-branes that describe AdS space 
sets the mass of the quarks. We consider the case where the bulk solution contains a black 
hole at temperature T. Although the probe branes are attracted to the black hole, at high 
enough quark mass the tension of the branes is sufficient to balance this attractive force and 
the probe branes lie outside the horizon in what is often called a 'Minkowski' embedding 



1 24, |2^]. In this phase, the gauge theory contains stable mesons with a discrete mass 
spectrum. However, above a critical temperature Tf, the gravitational force overcomes 
the tension and the probe branes fall through the horizon. The associated embeddings 
are referred to as 'black hole' embeddings. In this phase, the gauge theory has no stable 
mesons; the mesons melt. A detailed thermodynamic analysis of these embeddings reveals 
that the system undergoes a first order phase transition at T = Tf, where the probe branes 
jump discontinuously from a Minkowski to a black hole embedding pi, 25]. 



This phase transition has been studied extensively in systems of Nf probe D7-branes 
in the background of Nc black D3-branes, with Nf <C Nc- Without the D7-branes, the 
system is conformal and does not confine. Nevertheless, it can be a good model of a 
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confining system at sufficiently high temperatures T ^ T^. We now specialize to this 
model and fix the notation to be used below. 

The background metric for the black D3-branes may be written 

ds' = (If i-fdt' + dx') + ' + u'dnl^ , (4.1) 

where / = 1 — (uq/u)'^. The gauge theory directions are x = (x^,x^,x^), and the line 
element of the is 

dnl = d9^ + sin^ 9dnl + cos^ edcj)"^ . (4.2) 
The length scale L can be expressed as 

= At^QsN^II = 2Ui , (4.3) 

where Qs is the string coupling, ig is the string length, and we have defined the 't Hooft 
coupling A = 2TTgsNc. The horizon is at ti = uq, which is related to the temperature T 
through 

T=^. (4.4) 
As shown in |24, 25 1, the critical temperature for the meson melting phase transition is 

r;~M~^, (4.5) 

where M characterizes the meson mass gap at temperatures well below the phase transition, 
the quark mass is Mq ~ uq/I'^, and factors of order one have been neglected. 

Now, in the near horizon limit, the metric ( |4.1| ) can be written in the standard Rindler 
form by defining the coordinates 

ttTz^ = u — uo , y = L9, (4.6) 

and rescaling x ^x. Then in the limit of small z and small y, the background metric 
becomes 

ds'^ = -{2^Tfz^df + dz^ + dx^ + dy"^ + y'^d^l + dxl , (4.7) 
where x = (x^,x^,x^). For further details we refer the reader to |24, 25, p^ . 



4.2 Bubbles of melted mesons 

Consider now a thermal state of the gauge theory at some temperature T near Tf but 
satisfying Tf > T. The probe branes in the gravity dual lie on a Minkowski embedding 
which we may characterize by the minimum height zq of the branes above the black hole. 
This minimum height occurs at y = and is related to the constituent and bare quark 
masses M^Mq through [25| 

M, ~ ^ ~ iO-M„ (4 8) 

The constituent mass Mc is the physical quark mass taking into account thermal correc- 
tions. 
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From |25| we see that Minkowski embeddings with zq < 0.15 L are thermodynamically 
unstable and wih transition to black hole embeddings. Thus zq ^ 0.15 L. However, we 
wish to use the Rindler-space approximation ( [4.7D , which neglects terms in the metric of 
order 



u — Uq / z n 2 



Uq \L 

Thus for z ~ 0.15L, the error is of order (0.15)^. This occurs when T is close to Tf. While 
not parametrically small, the above error is certainly not large and we expect that the 
Rindler approximation should correctly capture much of the physics. Below, we restrict 
attention to such cases (which implies Mc ~ (0.15)^Mg) and use the Rindler approximation 
without further comment. 

Suppose that we now perturb this equilibrium by adding a massive quasi-particle which 
would have been stable (or at least very long lived) at T = 0. In the presence of the thermal 
bath all such quasi-particles become unstable and decay, but at least at large A, Nc some 
massive quasi-particles have lifetimes At much greater than their inverse masses m~^. 
In the gravity dual such resonances merely correspond to raising some object (graviton, 
massive string state, or black hole) to a height z ^ zq above the black hole horizon. We 
take the wave-function of this excitation in the gauge theory directions to be a Gaussian 
wave-packet centered on zero momentum with small width in momentum space, so that 
any spreading of the wave-packet is very slow. Due to the object's large mass, such wave- 
packets may nevertheless also be well-localized in position space. 

Since our object was raised to a large height z, it quickly begins to fall. By the time it 
reaches z ^ zo, the object is moving very rapidly. We restrict attention to an object which 
remains well- localized in the gauge theory directions (i.e., a beam- like graviton wave-packet 
directed downward, a massive string, or a small black hole), so that we may approximate 
its effect on the probe brane by that of an Aichelburg-Sexl shock wave of the form discussed 
in section |^ above. 

Now, there are two settings that one might consider, depending on whether the object 
is localized or delocalized on the S^. In the particular model considered here, the 
has size of order L. Due to various dynamical instabilities such as the Gregory-Laflamme 
instability for black holes one expects stable objects that are well- localized in the 
gauge theory directions and which satisfy r ^ L to also be well- localized on the S^. 

However, this fact means that the wavefunction of such objects will have components 
with significant charge under the SO (6) i?-symmetry associated with this sphere. If one 
eventually wishes to model more QCD-like systems without an iZ-symmetry, then one may 
expect that the case without i?-charge provides more universal results. Thus, despite the 
fact that such solutions are unstable in our model, we concentrate below on Aichelburg- 
Sexl-like solutions which are homogeneous over the S^. 

In our Rindler-space near- horizon limit, this is merely the spanned by y, O3, x^. 
One may think of the resulting shock- wave solution as being described by a (i = 5 Aichelburg- 
Sexl solution ( ^ni ) at each constant value of y,0,3,x^. In the notation of section |2| above, 
one may simply take the transverse radial coordinate to be = x ^ and d = 5. 
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We found previously that null geodesies undergo a discontinuous shift when they reach 
the plane of the shock. In fact, the same is true of any worldline subject to finite forces, 
as such forces add only a finite right-hand side to equation ( |2.5| ). Since the shift diverges 
at p = 0, it takes at least a portion of the D7-branes' worldvolume behind the black hole 
horizon. A localized bubble of the melted meson phase has been formed. 

We wish to describe the size and dynamics of this bubble. A full solution would require 
solving the D7-brane equations of motion in a dynamic background, taking proper account 
of the boundary conditions at infinity. We can, however, make some progress within our 
Rindler approximation. Note first that, due to the scale invariance of Minkowski space, 
in our Rindler approximation the curvature scale of the initial static D7-brane at the tip 
{y = 0,z = zo) must be set by zo ||2^. Thus, for y <^ zo we can model the initial 
configuration of the brane by the hyperplane z = zq, with xg = 0. In terms of the light- 
cone Minkowski coordinates, the brane lies on the hyperbolic curve —UV = Zq with xg = 0. 

At V = Vq, a point on the brane at impact parameter p will be shifted by an amount 
|Af7| = fi^{p). Thus, immediately after the jump the profile of the branes is given by^ 




4p) = \I4 - iTZ^'^ip) = - — , (4.9) 

where 

Ro = -^ (4.10) 

vr Tz^ 

and G5 = Gio/('7r^L^). We will comment further on the fact that Rq ^ E in section |5|. 

Equation (|4.9| ) describes a cylinder of radius Rq and height zq which attaches smoothly 
to an infinite plane at z = zq. Our approximation is valid for Rq <C zq, in which case the 
cylinder is tall and narrow. It therefore has a tendency to contract. Furthermore, as shown 
in the appendix, over most of our Rindler region (for p < = {RqZq)^^^ /2, 01 z < z^ with 
z* « 2:0 — (-^o-^o)^''^) the bubble nucleates with dp/dV < 0. Thus, unless other influences 
from beyond can arrive in time to prevent it, this part of the tube will collapse. 

We now show that no such influences can arrive in time to prevent collapse. Some care 
is required due to the fact that our initial data is given on the surface V = Vo,z = z{p), 
which is not a constant time surface in any natural coordinate system. To begin, consider 
a slice of the brane at constant z. If no influences arrive to slow the collapse, this part of 
the brane will collapse at V = Vq + AV for 



*When the falling object is localized on the one obtains instead 



(R ^6 _ 4 Gio-E 
^ " 3^7 Tz^ ■ 



6 
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where in the last step we have used results from the appendix. Recalling that U = —z^jV 
we see that an event at which a point on the cylinder collapses to zero size has spacetime 
coordinates (V, C/, p) ~ (Vo + AV^, -zg/Vo + AVz^/Vo^, 0), where we have assumed ATZ/Vb < 
1. On the other hand, the event at ^ = Vq where dp/dV becomes positive has spacetime 
coordinates = Vo,C/* = — 2;^/K,P*) ~ (^o, — ^o/^Oi P*)- The spacetime interval between 
these events is therefore 

A,^ = -AJ/AV + (Ap)^~-(^-J^)%p^ (4.12) 

which is positive for /) < /? ~ -i^R^^ . Thus, the tube collapses at least out to height 



5/9 



It is interesting to ask about the rate at which the bubble collapses. However, this 
question is complicated by several considerations associated with the fact that the bubble 
is far from even local thermodynamic equilibrium. First, the bubble does not form at a 
single value of the gauge theory time t. Instead, it forms on the null surface V = Vq. 
Second, there is some ambiguity as to what one wishes to call the "size" of the bubble. 
The naive bubble size is Rq, the radius of the cylinder of D7-brane intersecting the horizon 
at V = Vq. However, we have seen that the rather larger region of the brane out to 
p = p ^ Rq eventually collapses. A causality calculation similar to the one above then 
shows that any excitation of the brane (e.g., a meson) which begins at p < p ~ 2(i?^zo)^/^ 
is necessarily caught in the collapsing bubble, so that it is the scale p which governs the 
region in which meson propagation is highly disrupted. In terms of the gravity side of 
the gauge/gravity correspondence, we will tentatively describe these effects as a bubble of 
proper size Rq surrounded by a rather larger region in which propagation of disturbances 
on the D7-brane is highly disrupted (though we again emphasize that in reality the entire 
region is far from any local thermodynamic equilibrium). Using the metric (4.1) to translate 



distances to the gauge theory, for small Rq/zq we have a bubble of size Rq = Rq/{-kTL) ~ 
E / {7r{0.15)'^ N^T^), surrounded by a large halo that interferes with meson propagation. 

On the other hand, the time interval At between the creation and collapse of the 
bubble is a well-defined quantity. Using V = ze^'^^* and AV/Vq <C 1, we find 

Ay 



2ttTVo 

^ f^y^'- 1 f ^^y^' forp=-p. (4.13) 



At 



2ttTzo AttT \ zo J 4(7rr)2/3 I 0.15 J 

When the tube collapses, the motion can no longer be described using the classical equations 
of motion for the D7-brane. However, one expects that the collapsed tube breaks off from 
the rest of the D7-brane, perhaps annihilating in a shower of radiation. What remains is a 
D7-brane still somewhat distorted from its equilibrium configuration, which then oscillates 
in z until the remaining energy is dispersed. 

Let us briefly comment on larger bubbles, beyond what can be accurately described 
in our Rindler approximation. Recall that zq > 0.15L so that Rq ^ zq implies Rq ^ L. 
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Thus, the dynamics of such bubbles will proceed slowly enough to be affected by the shape 
of the branes in the asymptotic region of the spacetime. It is just this asymptotic region 
which controls the difference in energy between the most stable black hole and Minkowski 
embeddings found in |24, p5[. Since we work at a temperature where the Minkowski 



embedding is thermodynamically stable, it is clear that the energetics of this region cause 
the bubble to contract; i.e., the bubble is compressed by thermodynamic pressure. 

While our Rindler approximation does not allow us to calculate the precise timescale 
for this collapse, the collapse is given by the motion of a probe D7-brane in the black 
D3-brane background. As a result, the collapse timescale must be of the form times a 
function of Rq, the initial height zq, and the AdS scale L. The multiplicative factor of T^^ 
comes from the metric ( [4.7D and appears when translating proper time to coordinate time. 
The proper time, however, cannot depend on the temperature T since, in the particular 
model considered, one may change T via a change of coordinates. In gauge theory terms, 
this is associated with the scale invariance of the = 4 super- Yang- Mills theory to which 
the quarks have been added as a small perturbation. From the gauge theory perspective 
the timescale is determined by T, Tf , and the bubble size Rq and is otherwise independent 
of X,Nc,Nf. Because additional dependence on T appears when expressing zq in terms of 
gauge theory parameters. At need not be proportional to T~^. 

It is also interesting to analyze the bubble dynamics directly in gauge theory terms. 
In contrast to the above, we now assume that the 't Hooft coupling A is small in order that 
we may calculate reliably. 

In the gauge theory, the bubble of melted mesons is distinguished from its surrounding 
plasma by the presence of deconfined quarks. The bubble shrinks due to processes where 
these quarks annihilate to form gluons which may then leave the surface of the bubble 
and carry away energy. Quarks can also combine to form mesons, but this process is 
suppressed by a factor of 1/Nc in the large Nc limit. The quark-quark-gluon vertex factor 
is gvM = V^A/iVc so the production rate scales as 




^f^cixhr) (4.14) 



in the limit of large Nc- Here the factor NjN^ comes from summing over initial states. 
Since there are Nj species of heavy quarks, and since each species comes in Nc colors, the 
lifetime of the bubble scales as NjN^. This result matches the scaling with Nc,Nf found 
on the gravity side at large A, where collapse of the bubble was merely classical brane 
dynamics in a fixed background. 

One may expect the bubble to decay only due to interactions near the surface so that 
the lifetime of the bubble scales with the volume-to-surface area ratio Rq. This expectation 
is justified when the interior of the bubble is in local thermal equilibrium so that it is equally 
likely there for two quarks to combine into a gluon, or a gluon to decay into two quarks. 
However, since it exists only for a time At <^ ^/T, our bubble is in fact quite far from 
thermal equilibrium and it is possible that the entire bubble contributes to the decay of 
the bubble. 
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5. Discussion 



The above sections used gauge/gravity duality to describe two types of processes involving 
the rapid transfer of localized energy to deconfined plasmas at large \,Nc- The first 
process involves high energy collisions with balls of deconfined plasma surrounded by a 
confining phase, while the second involves the localized heating of a deconfined plasma. 
Both processes correspond to collisions with black holes in the gravity dual, where they 
result in the nucleation of a new "arm" of the horizon that reaches out in the direction of 
the incident object. In the first case, the arm stretches into a gauge theory direction while 
in the second it extends in the holographic direction. We have described the shape of this 
arm in the limit where the incident object is highly boosted and where one considers only 
the region close to the original black hole. After the collision, the arm is absorbed into the 
body of the black hole, though we have not described this process in any detail. 

Within the context of this duality, our results show that highly boosted particles 
incident on a plasma ball induce the formation of "virtual arms" which, though they 
contain no plasma themselves, will with large probability cause any additional particles 
passing through them to be captured by the plasma (see Fig. |^. These virtual arms 
summarize the interactions of both the original plasma ball and the incident highly boosted 
particle with any additional particles that may be present. The shape of a virtual arm was 
determined in section ^ (see equation (^]^)), where the arm length was found to scale 
with the energy of the incident particle and inversely with a power of the cut-off f on the 
transverse size. 

It is straightforward to generalize our single-particle analysis to the case of a dense 
cluster of particles, such as might be present in a pulsed beam. There are two essential 
points: i) The arm shape is driven by the Aichelburg-Sexl potential which is determined 
by solving a Poisson equation whose source is the incident energy. Thus, this potential 
satisfies a superposition principle, ii) Since the dynamical timescale for the horizon is 1/T, 
the response of the horizon is independent of the distribution of the incident energy along 
the direction of motion as long as the entire pulse of energy arrives in a time At <^ 1/T. 
Thus, (3^) may be used to calculate the length of the virtual arm induced collectively by 
particles in the cluster by taking the energy E to include all incident particles that arrive 
within the relaxation timescale 1/T and r to describe the size of the incoming cluster of 
particles. In fact, if the cluster is rotationally invariant about the direction of motion, ( |3.5| ) 
can be used to compute the arm shape by taking E(r) to be the energy contained within 
a cylinder of radius r. 

The scaling with E is universal in the limit of a thin relativistic beam, independent 
of the detailed model for the duality. The scaling with the cut-off r depends only on the 
dimension of the gravity dual. It would therefore be interesting to compare both scalings 
with calculations in a weakly coupled gauge theory and with experiments that might one 
day be performed with QCD plasmas. Such experiments could provide a new test of the 
extent to which QCD plasmas might be described by a classical gravity dual, independent 
of the details of any particular model. 

We also obtained results in the context of local heating of plasmas. Not surprisingly, 
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we found that such heating can produce bubbles of a higher temperature phase, but that 
these bubbles are unstable and collapse. In terms of the proper bubble size Rq in the bulk 
or the proper bubble size Rq in the gauge theory, we were able to study the dynamics in 
detail for T ^ Tj, large A, and small bubbles Rq <^ zq, Rq <ti 1/T. There we found that 
the bubbles are surrounded by a much larger region where meson propagation is highly 

1 / B 

disrupted. We also found that the bubbles collapse within a time At ~ ^^^^^2/3 ( cnk ) 
More generally (for larger Rq, smaller T, or small A) the collapse time may depend on 
T, Tf, and Rq, but must remain otherwise independent of A, Nj, Nc- 

A result of particular interest is the dependence of the bubble radius on both the 
energy E and the temperature T. Recall that we considered an example four-dimensional 
gauge theory with Nf flavors of fundamental matter in which the mesons were stable at 
T < Tf, but where the bubble corresponded to a phase in which mesons are melted. It is 
clear that little of the incident energy is transferred to the mesons, so the main effect is a 
local heating of the plasma. Since the specific heat of the system is of order A'^^ and the 
mesons melt at temperature Tf, one might therefore expect the bubble size to be of order 

Rnaive ~ I (^rp^ _ X)T^N'^ j ^^''^^ 

as there is only enough energy to raise the temperature from T to Tf over a region of this 
size. 

However, at least when T Tf, the bubbles that form in our processes are in fact of 



size Rq ~ 7i-(o i5pT^Af2 ' where we have translated ( 4.10|) to gauge theory parameters using 



Rq = Rq/{ttTL), zq ^ L and ([4. 3D. Thus, the bubble size is independent of T — Tf and, 
in the limit of large E, the actual bubble size Rq is much larger than Rnaive- This is due 
to the fact that the bubble is quite far from thermal equilibrium. In fact, from the bulk 
perspective the black hole arm that pulls the probe brane behind the horizon exists only 
for the time scale At ~ 1/T typical of thermal fluctuations. Thus, there is no reason to 
assign the arm an energy density characteristic of thermal equilibrium at temperature Tf. 

Our analysis considered only physics associated with the region close to the original 
black hole horizon and cases where the incident particles approach the black hole at rel- 
ativistic speeds. In such regimes, the behavior of the black hole horizon is universal and 
independent of the details of any particular duality. One may therefore hope that this limit 
describes universal physics of deconfined plasmas at large Nc- However, it is also clear that 
there is much to learn from the details of particular models in regimes where these limits do 
not apply. Some such extensions of this work may be straightforward, such as taking into 
account the details of the black D3-brane solution and the various D7-brane embeddings 



found in [24, ^] . Other more complicated extensions might involve studying the collisions 
of two plasma balls. Such situations are also more complicated from the bulk viewpoint 
and may require sophisticated numerical simulations. However, recent progress (see e.g., 
||2^ ) in simulating black hole collisions in 3+1 dimensions suggests that such calculations 
may nevertheless prove tractable in the near future. 
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A. Brane dynamics in the Aichelburg-Sexl solution 

In this appendix we provide a brief calculation of the D7-brane's initial velocity induced 
by crossing the Aichelburg-Sexl shock wave. We find that there is a region between z = 
and z ~ zq ~ (^o-^o)^''^ where the bubble nucleates in a contracting phase, as claimed in 



section 4.2 



In analogy to the solutions for null geodesies found in section 2.1, let us assume that 
the D7-brane configuration is described by specifying U(y,p). Hence, the induced metric 
on the brane takes the form 

-— + pHp)5iV - Vo)] dV^ + -Q^dpdV + dp" + p'd^l + dy^ + y^d^l . (A.l) 
As a result, the Lagrangian density appearing in the DBI action for the brane is 



'dU 1 [duV 



^"^V^^n^^J -/.<^(p)5(F-^o). (A.2) 



Based on the discussion in section we choose an ansatz for U(y,p) of the form 

U{V, p) = -y + HV, p)e{V - Vo) , (A.3) 

where h{V, p) is continuous in V and satisfies h(yo,p) = \ AU\ = p^{p). Using this ansatz 
and the equation of motion obtained from ( |A.2| ) , we find that immediately after the shock 
wave passes 

du _ zl dh _ zl p^^\pf 



dV ^ dV 



V=Vo "0 

Now, recall that the Rindler coordinate z is related to U, V through z = \/—UV, so the 
velocity in the z direction (at fixed p) is 

The radial velocity of the cylinder (at fixed z) is then given by 

dp^ ( dz\ f dz^^^ if RoZq\ 



dvj, \dv ) p \dpjy Vo V ; ' 

where we have used z{Vo,p) = zo\/l — Ro/p and Rq is defined in ( [4.10 ). Equation ( [A.6D 
implies that dp/dV is indeed negative for p < p* = - — ^ — , though it is positive for 
p > p*- 
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